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Abstract 

We review to what extent charge independence breaking {GIB) in isospin related 
reactions of the type A^A^ — > dir can or could be seen in existing data. In doing this 
we present fits to global and threshold cross sections including most recent data. 
Applying these we point out the probable impossibility to make model independent 
predictions even based solely on the different threshold energies and the well known 
pp d-K~^. A possible discrepancy is seen between the np — > dn^ and pp dn^ 
data, which may require invoking explicitly isospin symmetry breaking interactions. 
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1 Introduction 



It is well known that the isospin invariance as introduced originally by Heisen- 
berg is only an approximate symmetry even in the case of strong interactions. 
While in nuclei it is relatively strongly broken, it is, however, very good on the 
level of nucleons, where the mirroring of neutrons and protons is related to 
charge symmetry (CS), a more accurate symmetry than charge independence 
{CI). In the framework of the meson theory of nuclear forces a classification 
of different isospin structures is given in Ref. [1]. After Weinberg's ground 
breaking work on the quark masses [2] it became clear that in addition to the 
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Coulomb force, the difference between the hght quark masses breaks isospin 
symmetry. This means that in QCD under exchange of the u and d quarks 
there will be physically observable changes. Therefore, measurements of the 
violation of the isospin symmetry might allow the deduction of this mass differ- 
ence. Recent reviews related mainly to this charge symmetry breaking ( CSB) 
are given in Refs. [3,4]. For a recent general comprehensive review on meson 
production we refer to [5] and [6]. 

Isospin symmetry breaking effects are small and few dedicated experiments 
have been performed to study them in inelasticities. Nevertheless, in this ar- 
ticle we attempt to draw together what is known about the relation between 
the two reactions 

np TT^d (1) 

and 

pp TT+d (2) 

paying special attention to their differences. It is readily clear that these are 
not related by the charge symmetry operation (a reflection of the isospin in the 
2; or direction), but by charge independence, a wider symmetry of arbitrary 
rotations in the isospin space. 

As pointed out by Yang [7] already in 1952 isotopic spin invariance requires 
that the cross sections must be related by 



daipp — > TT'^d) 
da{np — > TT^d) 



2. (3) 



This can be shown by expressing the initial and final states in terms of eigen- 
states of the total isospin (spin summation implied) 

dajpp ^ TT+d) ^ 

dainp^Ti'^d) |1/V2 (1, 0|5|1, 0) + 1/^2 (1, 0|5|0, 0)|2 ^ ^ 

where S is the scattering matrix. If the isospin is conserved as required by 
charge symmetry, the second term in the denominator must vanish. Also, if 
the strong forces are independent of the third component of isospin, i.e. charge 
independent, then 

(1,1|5|1,1) = (1,0|5|1,0), (5) 
and the amplitudes in the numerator and denominator cancel leaving the 
factor of two. 

CSB requires typically a vector operator (tensor of rank one) in the isospin 
space as explicitly seen in the classification of two nucleon forces in Ref. [1], 
i.e. the isospin structures (ti + T2)o (class 111) or (ti — T2)o or (ti x T2)o (class 
IV). Here the latter operators necessarily change the isospin and thus also in 
the np — > dTT^ reactions (and only there) can mix isospin zero initial state 
in the production amplitudes. This causes an asymmetry of the cross section 
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about 90° in the centre of mass system. Otherwise the cross section would be 
symmetric [8] . A measurement of the charge symmetry breaking asymmetry in 
the meson production reaction np — > 7r°d close to threshold has been published 
recently [9]. 

Charge independence is a broader symmetry encompassing any combination 
of charges, i.e. a symmetry with respect to arbitrary rotations in the isospin 
space. Isotensor class 11 forces of Ref. [1] (3rior2o — ti ■ T2), while respecting 
charge symmetry, do violate charge independence in the two nucleon case, 
and similarly again Eq. (5) need not be valid in meson production either. A 
deviation from the factor of two would be a sign and a measure of violation 
of isotopic invariance. It should be noted that in the present case and in 
NN scattering the isotensor operator cannot change the isospin, which has 
to be unity in both the initial and final state. Therefore, the spin-spatial 
symmetries would be the same as in the isospin symmetric case, contrary to 
charge symmetry breaking. 

Unfortunately even less is available for charge independence breaking (CIB) 
than CSB in meson production, although, in principle, this should be signifi- 
cantly larger than CSB, a few % vs. a few thousandths. To our knowledge, only 
the unsuccessful attempt [10] to see CSB reports possible marginal CIB in the 
angular dependence at 795 MeV. The experimental difficulties here stem from 
two origins: Firstly, there is not a "null experiment" nor such a clear forward- 
backward asymmetry signal as in CSB. Secondly, these measurements require 
clearly a comparison of separate experiments, one with neutron beams, prone 
to meet normalisation problems. 

However, further motivation to study isospin symmetry breaking for the above 
reactions (1) and (2) is given by the recently measured, somewhat related more 
complicated reactions 

pd 7r° ^He (6) 

and 

pd ^ 7r+ -^H . (7) 

These reactions have the advantage that they can be measured simultaneously 
and thus their ratio does not suffer from normalisation uncertainties. Indicative 
of isospin symmetry breaking, the excitation function of the ratio between 
the differential cross sections for maximal momentum transfer to the pion 
showed deviation from the isospin ratio two when passing the pd ?7^He 
threshold [11]. Full angular distributions for both reactions at lower energies 
were reported in Refs. [12] and [13]. There the angular distributions also show 
a slight but systematic difference. However, the origin of this difference could 
be different nuclear wave functions (possibly due to the Coulomb force) of ^H 
and ^He rather than CIB production mechanisms. It would be interesting to 
see if and how isospin symmetry breaking would persist in the more basic and 
theoretically cleaner reactions (1) and (2) without this distortion. 
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Fig. 1. Relative differences between the Gamow factors (see Ref. [16]). 

The basic theoretical arguments for handhng the data are presented in the 
next Section and then the results of the analyses in Section 3. 



2 Coulomb corrections and kinematic considerations 



Differential as well as total cross sections of the two reactions obviously differ 
because of the mass differences of the nucleons in the entrance channel and of 
the pions in the exit channel, and also due to the presence of the Coulomb force 
in the reaction (7). Furthermore, charge dependent forces may be involved in 
the reactions. Typically these would arise from charged vs. neutral meson mass 
differences in meson exchanges. Charge asymmetry, in turn, can get contribu- 
tions from nucleon mass differences, meson mixing [14] , magnetic interactions 
and explicitly isospin symmetry violating pion- nucleon rescattering [15]. 

In comparisons of the reactions (1) and (2) we shall first study the second 
effect. A common approach to correct for the Coulomb effect is to apply the 
s-wave Gamow penetration factor 

l^o(OI =^5^ with e = (8) 

to the cross sections. Here q is the centre-of-mass momentum of the pion and 

J] = q/mT^+. ^ 



^ One should note that the final state dependence of the partial wave amplitudes 

is through the pion momentum q and that in r/ the charged pion mass acts only as 
a unit and should be kept the same in both reactions (1) and (2). Use of different 
scaling masses would cause an artificial effect based on the choice of units. Therefore, 
the use of "natural" mass 171^^0 in reaction (1) is an unfortunate choice. 
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Fig. 2. Kinematical limits of the two reactions, tmax and tmin are shown as function 
of s. 

Niskanen and Vestama [16] have calculated more generally such factors under 
various assumptions: for a point-like charge distribution with pion angular 
momentum / (|Cip) and for an extended charge distribution (|Cf p) applying 
the Reid potential [17] to obtain the deuteron wave function. It was seen 
that the extended and point sources gave very similar corrections so that 
the use of an extended source is not really warranted with the experimental 
accuracy expected in near future. In Fig. 1 we show the relative difference 
of the Gamow factors for p- waves with extended charge |CJ^P and s- waves 
|CoP relative to |CoP as a function of the dimensionless centre of mass pion 
momentum r] = p-„/mTj-+. The first difference is large only in a range for very 
small momenta rj < 0.05, where p-wave cross section is negligibly small, and 
then the difference between the two factors decreases below 2%. This is the 
order of magnitude of precision that can possibly be obtained. 

Justified by these results in removing the Coulomb, for the sake of simplicity, 
we will make use of only |Co|^ jC^p 1.008 |Cop, since a separation into 
different partial waves is far from straightforward [16, 18]. As can be seen from 
Fig. 1 in relative angular dependence this average factor cancels off to a good 
approximation for extended sources, since |CiP and |CoP are within 0.5% 
from each other. The similar factor |CoP is also applied for the pp initial state 
giving a normahsation effect of 1-2%. 

We proceed now to study the effect of kinematical differences due to differ- 
ent masses. The differential cross section da/dt depends on the Mandelstam 
variables s and t. Here we use t as the four momentum transfer squared from 
the projectile to the pion. In Fig. 2 the kinematical limits, i. e. maximal and 
minimal values of t as functions of s are given for the two reactions (1) and 
(2). 

In addition to the different threshold energies ^/so also for given CMS excess 
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Fig. 3. Kinematical limits as function of the excess energy e = s/so- The t-scale 
is shifted for reaction (1) by —0.01065 GeV^. 



0.64 




0.0 0.1 0.2 0.3 0.4 0.5 

1 



Fig. 4. Same as Fig. 3 but as function of the dimensionless pion momentum 



energies e = ^/s — ^/so the t values are shifted to a good approximation by a 
constant 0.01065 GeV^ as can be seen from Fig. 3, where the hmits for the 
reaction (1) have been brought on top of the reaction (2). As a function of the 
final state momentum the rounded bottom of the lopsided parabola gets a V 
shape in Fig. 4. 

One may well argue that a large part of the difference between the two reac- 
tions comes, in addition to the Coulomb effects already described, simply from 
the different thresholds. This means that for a given final state momentum 
(the decisive dependence close to thresholds) the incident kinetic energy is 
slightly different. This gives a trivial effect in the phase space factor P (with 
an obvious notation for the energies and momenta p* in the centre-of-mass 
system, N — p or n, and including also the spin factor 1/4) 
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present in cross sections, but it influences also the transition matrices in R 
themselves. In particular, due to the lower threshold as a function of the final 
state momentum the amplitudes for np — > Tr^d should lag behind pp — > ir'^d 
in reaching and passing e.g. the A resonance. (One should note that the mass 
difference of the nucleons does not directly affect the proximity of the A, since 
similar mass differences would be present in that, too.) 

In Ref. [16] changes in the total and differential cross section were expressed 
directly by differences in both the phase space and spin amplitudes as 

Sa ^ 2a{np) - a{pp) = ^"^ + ^"" (i?., - Rpp) + (P„, - Ppp) ^"^^^^^ . (11) 

Results for the relative (integrated) differences in P and R were given, for the 
latter by a model calculation. 

One might choose another path aiming at a model independent prediction for 
the difference based on existing data. First dividing the trivial phase space 
factor off directly from the empirical data as well as the Gamow factors (8) 
one can study the energy dependence of the more interesting R for equal pion 
momenta in the two reactions. Then one would get from 

H R 

SR = Rr^p -Rpp^-^ {E,{np) - E,{pp)) (12) 

a prediction for the dynamical difference between the reactions (1) and (2) 
from data on only the latter. However, if the quantity Rpp is given as a function 
of 7] like usual for the cross section, a singular drj/dEi arises. Consequently 
this product close to threshold is too large for the first order estimate to be 
useful. 

It would now be tempting to think that deviations from this prediction would, 
in principle, be a direct indication of charge dependent forces. However, this 

expectation is based on the kinematic relation between Ei and 7] used to 
connect the two different reactions. It still leaves untouched the possible dy- 
namical difference in the initial states due to the change of the energy but 
without any change in the final momentum. This will be seen later in Fig. 7. 
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Fig. 5. Comparison of differential cross sections for reaction(l) shown as squares 
and for reaction (2) shown as dots. The sohd and long dashed curves are param- 
eterisations discussed in the text. The short dashed curve is the SAID-fit [21] at 
r? = 0.51. 

3 Data Analysis 



We are now ready to compare data for the two reactions of interest. Unfor- 
tunately there are no differential cross sections for the same 77 or e and the 
same cm emission angle. Nevertheless, in order to stress this shortage we com- 
pare two angular distributions, at reasonably similar energies. Fig. 5 shows 
the centre-of-mass differential cross sections from the data of Wilson et al. for 
the reaction (1) [19] together with those for the reaction (2) from the GEM 
collaboration at COSY [20] (reflected about cos ^ = 0) . 

The differential cross sections can be expanded in terms of Legendre polyno- 
mials P;(cos^) [22] 



47r^^^%^ = ^(7?)Po(cose) +^(7?)P2(cos^) (13) 
all 

^airj)[l + a{rj)P2{cose)], (14) 
with the indicated normalisation 

a{r)) = Aoiv) (15) 

and the anisotropy 

a{r])=A2iv)/Aoiv)- (16) 

Eq. (13) is valid for both the reaction pp n^d and for reaction np Ti^d 
with the same coefficients, if charge independence is obeyed (apart from the 
normalisation factor 2). As an example this functional dependence is fitted 
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Table 1 

Uncertainty of the fitted parameters Eq. (13) to the data shown in Fig. 5. 



reaction 


data 


1] 


AAo/Ao 


AA2/A2 


Aa/a 


np — ^ 7r°c? 


[19] 


0.59 


0.027 


0.064 


0.048 


PP — > TT'^d 


[20] 


0.51 


0.0076 


0.023 


0.018 



to the data in Fig. 5. The uncertainties of the fit parameters are compiled in 
Table 1. 

The error in the charged total cross section would be small enough to see 
isospin breaking effects, while that for the np reaction is of the same order 
as those effects computed in Ref. [16]. Furthermore, one should note that 
the np data are not absolute in their normalisation. The error in A2 is still 
significantly larger. However, interestingly the relative error of the anisotropy 
in the form (14) is smaller than what would be obtained from the form (13). 
Apparently this reflects the fact that normalisation uncertainties present in 
A2{i]) are not present in the relative 0,(1]) making the latter a preferred form 
to fit. Of course, the actually obtained fitting result is the same in both. The 
anisotropy will be discussed in more detail further down. It is also of interest 
to note that the SAID phase shift analysis [21], which does not include the 
data from Ref. [20], is significantly different. 

In the comparison (Fig. 5 and Table 1) we have not applied the corrections 
discussed above, because of the differences between the tw^o data sets. It is 
evident that the precision and quality of cross sections in proton-proton in- 
duced reactions is far superior to those from neutron-proton induced reactions, 
where high quality data are still lacking. Therefore, in the following consid- 
erations we do not use individual data sets but rely rather on more global 
parameterisations of almost all available data as functions of the pion angle 
and momentum. 

We now consider the total cross section more globally as function of 77. In Fig. 
6 we show almost the world data set of cross sections for pp —>■ ir'^d (Refs. 
[20,23-41]) After a steep rise the cross section falls down after the maximum 
which corresponds to the centre of the A resonance. A global fit performed by 
Ritchie [42] underestimates the new near threshold data. Therefore we have 
performed a new fit to the data body after removing the phase space and 
the Coulomb effects. To get a reasonable fit to the partially contradictory 
data it was found necessary to make some choices concerning their inclusion. 
The following data were omitted in the fit: The data from Ref. [43] (see the 
discussion below). Two points from Ref. [28], one point from Ref. [44] in the 
interval 0.661 < rj < 0.735, and the lowest point from Ref. [27]. 
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Fig. 6. Upper part: The total cross section for pp 7r~^d as function of the dimen- 
sionless pion centre of mass momentum. Lower part: Same as the upper part but for 
the angle integrated sum of the squared matrix elements, Rpp (Coulomb corrected). 
The solid curves are fits which are discussed in the text. 
Table 2 

Fitted parameters of Eq. (17) to the total cross sections. 

a(mb) b c d e 

0.17 ± 0.03 1.34 ± 0.06 1.77 ± 0.04 0.38 ± 0.02 0.096 ± 0.02 



We find the function 



arj 



(c - e*0 + e ^ ^ 

to account for the raw cross section data. This fit is included in the upper 
panel of Fig. 6 and the fit parametres are given in Table 2. 

Also shown in Fig. 6 are the Coulomb corrected angle integrated matrix ele- 
ments. These can be excellently fitted by the squared Lorentz function 



R 



■pp 



1 2 



(&2 



(18) 



where 62 sets the resonance position and 63 the half-width. The corresponding 
curve is included in the figure and the fit parameters are given in Table 3. We 
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Table 3 

Fitted parameters of Eq. (18) to the angle integrated and Coulomb corrected matrix 
elements. 

61 (fm^^MeV) 62 63 

1252.347 ± 11.232 1.2503 ± 0.0064 0.7922 ± 0.0108 
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Fig. 7. Predictions for the relative difference of the total cross sections of the re- 
actions (1) and (2). The dashed curve is the model prediction of Ref. ( [16]) for 
SR/Rav, while the solid curve is based on Eq. (12) with Rpp . 

also tried adding a constant and ry-dependent background, but the quality of 
the fit remained essentially the same and the additional terms were consistent 
with zero. 

Although the above form (especially with background terms) may look super- 
ficially like a theoretically motivated squared amplitude, this is not in reality 
the case. Since different spin amplitudes do not mix in the cross section, odd 
pion waves (starting with the p waves) should have a separate squared am- 
plitude with overall threshold behaviour rj'^ (and the A(1232) resonance) and 
even waves their own with a constant threshold behaviour (to which even 
powers of r] would be added at the amplitude level). The above form is only 
justified as an attempt to fit the resonant structure with a simple function and 
it is, indeed, remarkable that nearly all total cross section data can be fitted 
with just three parameters. One may note that this form is not even analyt- 
ically of the correct form, which should be a function of 1]^. With functions 
depending on only this we were not able to get fits of comparable quality even 
by increasing the number of parameters. The data appear to require a linear 
term in 77 in particular in the neighbourhood of ~ 1. Also theoretical treat- 
ments have had trouble in giving enough cross section in this region without 
destroying agreement elsewhere. With the threshold behaviour given by this 
form (i.e. dR/drj as — > 0) the prediction (12) becomes excessively large 
as can be seen from Fig. 7. 
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In Eq. (12) the above functional form is used to make a prediction for the 
difference between the reactions (1) and (2). This is shown together with the 
model prediction of Ref. [16] in Fig. 12. The agreement between the two is 
quite reasonable at higher energies. Considering the discussion of the weak 
model dependence in that reference the difference is significant. However, it 
cannot be attributed to charge dependent interactions (both are independent 
of those). As discussed in the previous section, the present prediction is based 
on the relation between the initial energy and final pion momentum (albeit 
with different thresholds) employing the known pp data. Possible difference 
in the initial state at different energies does not come automatically except 
through the final momentum dependence. This difference is most pronounced 
for small values of rj. Even with the correct threshold dependence the uniformly 
increasing matrix elements should predict smaller np than pp cross sections. 
However, the converse is seen in the model prediction (dashed curve). Appar- 
ently, if one could have a varying initial nucleon energy with a constant final 
momentum, the matrix elements would decrease with the energy. This trend is 
overcome only in the region of the sharp rise of the A-resonance effect, where 
the np reaction lags behind and the cross section becomes smaller than that 
of the pp initiated reaction. 

We now compare the angle integrated matrix elements as function of 77 in the 
threshold region where there are rather precise data for both reactions. Here 
the cross sections for np — > ir^d are taken from Ref. [45]. They are in the 
threshold range up to 77 = 0.32. In this range one should be able to apply the 
low energy expansion 

ao(l + a;i7?2), (19) 

where the first term is predominantly s-wave and the second p-wave [18]. 
The Coulomb corrected cross sections for pp — > i^'^d are from Refs. [20, 23- 
27]. Here we have restricted the range to 77 < 0.5 to be well-matched. It 
is worth mentioning that the cross sections from Refs. [27] and [26] agree 
nicely with each other when the same Coulomb correction is applied. The 
Coulomb correction used in Ref. [27] is somewhat too large (see the discussion 
in Ref. [16]). Here we have started with the uncorrected cross sections from 
Ref. [46] and applied the corrections discussed above. The deduced matrix 
elements are shown in Fig. 8 and the corresponding fit parameters in Table 4. 
Those for np — > n'^d reactions are multiplied by the isospin factor of two. Fits 
with Eq. (19) are shown as error bands. 

The excitation functions for the two reactions clearly differ from each other 
with respect to magnitude and shape. Applying the same corrections as for the 
data we have also extracted the matrix elements from an earlier fit [42] to data, 
which appeared before the advent of new data in the threshold region [25-27]. 
This is also shown in the figure, where the different slope to the data shows the 
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Fig. 8. Fit of the low energy function (19) to the sum of squared matrix elements 
extracted from the data. The error band represents the confidence interval on the 
95% level. Data for reaction np —>■ n'^d, corrected by the isospin factor of 2, are 
shown as squares, those for reaction pp — Tr+d as dots. Also shown is the result for 
the fit from Ref. [42]. 

importance of the new threshold data. Following Feynman's last point rule^ 
in fitting Eq. (19) to the data we have excluded one point from Ref. [23] which 
is far from the band of the other data. The same is true for the points with 
smallest r) from Refs. [27] and [45] (see Fig. 8). 

Differing from other np — > n^d measurements, the normalisation performed 
by Hutcheon et al. [45] is to the simultaneously measured np pn scattering 
cross sections. The absolute value of this reaction was taken from the SAID 
phase shift analysis [48]. The other measurements normalized their counting 
rate to pp TT+rf cross sections. However, even the np pn data might not 
be completely independent of an assumed isospin symmetry (see Ref. [49] and 
Ref. [50]). If that would be so we would expect two bands running in parallel 
with only different normalisations. Unfortunately at the present time we are 
not aware of such np reaction data which do not suffer from normalisation 
problems to such an extent as to make comparisons meaningless. 

In spite of the above normalisation problems it should be noted that at thresh- 
old Rnp > Rpp as predicted by Ref. [16] (the dashed curve in Fig. 7). However, 
the experimental difference is larger and goes through zero much faster than 
the prediction. This difference of the only comparable data shown in Fig. 8 can 
be further elaborated into the form of the prediction of Ref. ( [16] in Fig. 9. 

^ According to Feynman [47] "there is a principle that a point on the edge of the 
range of the data-the last point — isn't very good, because if it was, they'd have 
another point further along" . 
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Tabic 4 

Parameters from fits of Eq. (19) to the matrix elements in the close to threshold 
area. 



reaction 


ao (fm^MeV^) 


oil 


np — >■ TT^d 


(185 ± 2) X 10^ 


3.97 ±0.22 


PP TT+d 


(144.5 ± 4.7) X 10^ 


8.57 ±0.52 



Quite clearly the disagreement between experiment and theory cannot be re- 
moved even by renormalising the np cross sections. In the figure this has been 
done by moderate factors 0.9 and 1.1. This discrepancy might possibly be an 
indication of an effect of charge dependent forces and apparently deserves a 
more dedicated study. 
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Fig. 9. Comparison of the relative change between the model results of Ref. [16] 
(dashed curve) and the fit of Eq. 19 (solid curves) directly and with ao{np) renor- 
malised by factors 0.9 and 1.1. 



To avoid problems with the normalisation we now study the relative anisotropy 
a{ri), instead. The nice feature of studying this quantity is that the phase 
space correction and the isospin factor cancel exactly. As seen below, to a 
good approximation this is true for the Coulomb correction as well. However, 
before discussing the data we should first have a more detailed look into the 
meaning of this observable. 

Due to parity conservation always 

L{np) = ± 1 , (20) 

so that even isospin conserving pion waves are associated with initial odd- 
triplets and vice versa odd pion waves with even-singlet states, while for isospin 
breaking ( CSB) amplitudes the converse is vahd ^ . With the interference of 

^ Ref. [51] quotes additionally a superfluous and incorrect condition J = L{np) ± 1 
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opposite parities the contribution from CSB isospin mixing is odd in cos^ and 
docs not appear in the expansion with even polynomials, if both forward and 
backward angular ranges are represented with equal weights. Charge depen- 
dence due to isotensor forces does not mix these parities, so for this study it 
is sufficient and consistent to consider only angular dependencies symmetric 
about 90°. 

The differential cross section can be expanded in terms of the lowest partial 
waves a J (here up to the p-waves) as 

da{ri, cos 9) 1 



dfl 4 



(|ao|^ + |ai|^ + lasl^) + (jaaj^ - 2v^3fJ (aoa^)) P2 (cos6 

(21) 

By comparison with Eq. (13) one identifies now l/4(|aoP + |aip + |a2nas the 
total cross section and 1/4 (^|a2p — 2\/23fJ (aofla)) ^2- Here the ^5*0 p am- 
plitude ao is significantly smaller than the isotropic s-wave ai or the dominant 
p-wave ^ p enhanced by the A excitation with increasing energy [52]. In 
the case of a vanishing amplitude oq, ^2(^7) would be exclusively of a single 
p-wave without phase shift dependent interference. Also for a nonvanishing 
oo, if the same Coulomb phase is associated with the both p-wave amplitudes 
ao and 02, it will cancel off at this level of expansion and the result would be 
formally the same for both reactions. 

In order to be more sensitive to the p-wave we have to extend the momentum 
range. However, then also higher Legendre-polynomials and higher partial 
waves have to be considered in fitting the angular distributions and possible 
interference effects involving different Coulomb phases can arise, first the d- 
waves with the s-wave pions. Their partial wave dependencies will not be 
discussed here. For the general relation between the Legendre-co efficients and 
higher partial wave amplitudes (also including many spin obscrvablcs) see 
Ref. [53]. Again both reactions should be well matched in 1]. In the upper 
end of energies lack of data from np n^d above 7] = 2 limits the range of 
comparison with pp — > 7T~^d. 

The data for reaction (2) arc from Refs. [20,24-27,31,34-36,54]. The data 
from [32] and [55] were excluded, since they are far above the other data. The 
data from the time reversed reaction from [43] were also excluded from the fit 
because of the large spread within this data set (see the discussion in Ref. [20]). 
In the case of reaction (1) only the data from Refs. [45,51,56] were taken into 
account. Those from [19] and [57] were excluded on similar grounds as in the 
case of the charged channel. Although the error bars are small (see also Table 
1) there seem to be systematical uncertainties. This is demonstrated in Fig. 
10. 



for initial triplets, which would prevent e.g. pion s-wave production. 
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Fig. 10. The ratio a = A2I for the reaction n-p — ^ 'K^d as function of 77 = pT^/m.,^. 
The sohd curve is a fit as discussed in the text. 
Table 5 

Parameters obtained from fitting Eq. (22) to the anisotropies. 



reaction 


a 


P 


7 


np TT^d 


1.114 ±0.007 


1.044 ± 0.017 


1.020 ±0.029 


PP — > Tr~^d 


1.111 ±0.007 


1.003 ± 0.009 


0.926 ± 0.020 



In order to make a comparison between the two reactions possible, we have 
fitted to data a function which describes the dependence very well, although 
it is not motivated by any theory. The function is 



a{ri) = a exp 




(22) 



The fitted parameters are compiled in Table 5. While the values for a agree 
with each other, the values for /3 and 7 are slightly different. This also reflects 
in slightly different momentum dependencies of the anisotropy as shown in 
Fig. 11. 

In order to study further the deviations of the anisotropy for the two reactions 
from each other, we plot the difference of the anisotropies in Fig. 12. Also the 

uncertainty is shown. We can state that a^p is larger than a„p in the vicinity 
of the threshold {rj < 0.5) and smaller for the range above, although the error 
bar is large. 



4 Discussion 



We have discussed the total and differential cross sections of the reactions 
np — > TT^d and pp — > 7r+d as functions of the pion momentum in order to allow 
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Fig. 12. The difference of the anisotropy for np ir^d and pp n'^d (thick curve). 
The confidence interval (on 95% level) is indicated by the shaded area. 

for a comparison with respect to isospin symmetry. This takes into account 
automatically the different thresholds but with the cost of having different 
initial kinetic energies. However, this difference can be dealt with easily to first 
order once Coulomb effects have been corrected for and the slight difference in 
the phase space. It is found that the differential cross section data at individual 
energies for the reaction np — > 7r°d do not have the quality to make such a 
comparison useful. 

Having first made a global fit of the total cross sections and the associated 
effective matrix elements for future use, we then compared a number of total 
cross sections in the threshold region. This limitation is enforced by the lack 
of higher energy absolute cross sections of the np — > 7r°(i reaction. The result 
is that very close to threshold 2a{np — > n^d) > a{pp n^d). contrary to a 
simple "model independent" expectation (12) but in agreement with a theo- 
retical anticipation [16]. It was not possible to get both the size and steepness 
of the difference as a function of the final pion momentum to agree with the 
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charge independent theory. 

The cross section in the interval studied is mainly due to the transition 
^Pi — > ^Sis, which dominates the outgoing s-wave with the partial wave am- 
plitude tti, and ^ ^Sip in the p-wavc amplitude 02- There is also the 
externally indistinguishable contribution from the deuteron D state, so the 
shorter notation ^'''^^Lj — is well justified. Apparently the threshold nor- 
malisation is dictated by the s-wave amplitude (which yields the size order 
of the two cross sections dependent on the initial energy), while the p-wave 
determines the slope and thus the zero crossing point of the difference. Simply 
renormalising the np cross section has opposite effects in the threshold differ- 
ence and zero crossing in a comparison of theory and experiment and is not a 
solution to their discrepancy. 

We then compared the relative anisotropy for both reactions, where phase 
space corrections are not necessary and also the Coulomb effects are minute. 
It is mainly due to the outgoing p-wave with the A enhanced transitions 
— > p. The comparison of the anisotropy for both reactions shows that 
a{np — > TT^d) > a(pp — > ir'^d) although the uncertainty is large. The reason for 
this is again the poor quality of the np data. As an example we mention the 
overlap region around rj ^ 0.35 between the two data sets [45] and [56], which 
differ strongly, as is shown in Fig. 10. The problem of energy resolution of a 
neutron beam and the absolute calibration of its intensity may be overcome 
by changing the incident channel incident channel at the cost of increasing 
the number of participating nucleons. A quasi-free reaction pd — > n^dps with 
a spectator proton ps was recently shown to be feasible in a storage ring with a 
cluster target [58]. Also the reaction dp — > ii^dps with a fast spectator proton 
may be considered. 
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